The effect of nonmagnetic impurities on the local density of states in s-wave 

superconductors 
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We study the effect of nonmagnetic impurities on tlie local density of states (LDOS) in s-wave 
superconductors. Tfie quasiclassical equations of superconductivity are solved selfconsistently to 
show how LDOS evolves with impurity concentration. The spatially averaged zero-energy LDOS 
is a linear function of magnetic induction in low fields, N(E = 0) = cB /Hc2, for all impurity 
concentration. The constant of proportionality c depends weakly on the electron mean free path. 
We present numerical data for differential conductance and spatial profile of zero-energy LDOS 
which can help in estimating the mean free path through the LDOS measurement. 

PACS numbers: 74.25.Jb, 74.25. Bt, 74.25.0p 



I. INTRODUCTION 

Since Hess and co-workers^ succeeded in measuring the 
LDOS in the superconducting NbSe2 there have been 
many reports and theoretical studies on the electronic 
structure of the superconductor in the mixed state. The 
novel experimental technique introduced, scanning tun- 
neling spectroscopy, enables one to measure differential 
conductivity (DC) a{r,V) at various positions r and 
bias-voltages V . DC is closely related to LDOS of the 
superconductor {ks — 1), 



dE 



N{r,E) 



AT o 
iVo cosh^ 



E + eV 
2T 



(1) 



where cttv is DC in the normal state, e is electron charge, 
7V(r, E) is LDOS at the position r and energy E rel- 
ative to the Fermi level, and A^o is DOS at the Fermi 
level in the normal state. Only in the limit of zero 
temperature T ^ DC and LDOS are proportional: 
a{r, V)/aN = N{r, \e\V)/No. At finite temperature, DC 
is actually thermally broadened LDOS. 

It is clear that at low temperatures DC should follow 
the spatial structure of LDOS. Two prominent features 
should be mentioned. DC measured at the vortex center 
revealed a peak at the Fermi level (zero-bias peak) that 
well exceeds a^- This indicates that vortex core can not 
be viewed as being "normal" at least in clean supercon- 
ductors. The zero-bias peak in DC originates from the 
zero energy peak of LDOS at the vortex center, which 
is due to the low lying bound states inside the vortex 
core. The other remarkable feature revealed in Ref. Q 
is a star-shaped DC around the vortex core measured 
at fixed bias-voltage, with star orientation depending on 
the bias-voltage value. The six-fold structure of DC in 
NbSe2 is coming either from the effect of the hexagonal 
vortex lattice, anisotropic s-wave pairing or anisotropic 
Fermi surface, and most probably it is coming from each 
effect simultaneously. Again, the star-shaped DC origi- 
nates from the star shaped LDOS in the vortex lattice. 



However, the measured DC does not follow the sharp 
features of the corresponding, theoretically calculated, 
LDOS even if the experiment is performed at very low 
temperaturci^ The height and width of the zero-bias peak 
was found to be sample dependent indicating impurities 
as a plausible explanation for the discrepancy. Indeed, 
impurities are inevitably present in superconducting sam- 
ples on which the experiments are performed. There- 
fore, it is important to quantitatively study how LDOS 
is changing with impurity concentration. This is the pur- 
pose of our paper. 

There is another one topic that we analyze in this pa- 
per: the effect of impurities on the specific heat field de- 
pendence. In s-wave superconductors low energy quasi- 
particles are trapped inside the vortex core. Therefore, 
zero-energy LDOS, spatially averaged, is proportional 
to the number of vortexes: N{E = 0) oc No^'^B, with 
B being the magnetic induction and ^ the size of the 
vortex core. This translates into the linear field depen- 
dence of low temperature specific heat given by Cg/T — 
2tt'^N{E = 0)/3. The nonlinearity in low field Cs{B) 
curve should be related to the gap anisotropy. In the 
case of anisotropic s-wave pairing, addition of nonmag- 
netic impurities can smear out the gap anisotropy, which 
can be tracked by examining Cs (B) curves. This kind of 
measurement has been performed on Nbi_3;Taa;Se2'^ and 
Y(Nii_2;Ptj;)2B2Ci^ The intention was to make the gap 
isotropic by adding impurities, but with the price to have 
rather dirty s-wave superconductor. The effect of impu- 
rities on LDOS notwithstanding it would be of value to 
study how field dependence of spatially averaged LDOS 
evolves with impurity concentration in the most simple 
case of s-wave superconductors. 

So far, the only systematic experimental study of the 
effect of disorder on LDOS is by Renner et alA In partic- 
ular they measured the zero-bias DC at the vortex center 
in the alloy system Nbi_j;Taa;Se2. Substitution of Nb by 
Ta leads to systematic decrease of the electron mean free 
path. On the other hand small changes in the electronic 
spectrum is expected since Nb and Ta are isoelectronic 
and with similar atomic radii. Zero-energy DC is found to 
be very sensitive to the impurity concentration. It grad- 
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ually disappears and for x = 0.2 the zero-energy LDOS 
in the vortex center is the same as that of normal phase 
iVg. It was even proposed that DC spectra can serve as 
a measure of quasiparticle scattering time. 

In this paper the problem of LDOS in presence of non- 
magnetic impurities will be studied within the quasiclas- 
sical equations of superconductivity. Quasiclassical ap- 
proximation is adequate in superconductors where coher- 
ence length ^ is much larger than the atomic length fc^^. 
LDOS is studied within the quasiclassical approximation 
by UUah et al.^ and Klein^ for the case of isolated vortex 
in the isotropic s-wave superconductor. The full self- 
consistent analysis of LDOS in the case of vortex lattice 
is performed by Ichioka et alfi' All these studies assume 
clean superconductor without impurities. Dirtiness of 
the superconductor is only roughly estimated by Kleini 
As for the single vortex case, the effect of impurities was 
studied in Refs. The effect of impurities on DOS in 

extremely high field, where the Landau level quantization 
of the electronic energies should be taken into account, 
is studied by Dukan and Tesanoviciii Those phenomena 
are beyond the scope of this text. As for the study of the 
extreme case, dirty superconductor, the reader is refered 
to Refs. [Tlia 

The paper is organized as follows. In section ^ the 
method of solving Eilenberger equations is presented. 
The readers not interested in technical details may skip 
that section. In section lllll spatial and energy dependence 
of LDOS and DC for various impurity concentration is 
shown. In section Hvl the effect of impurities on the spe- 
cific heat field dependence is discussed. 



II. METHOD OF SOLUTION 

There are various methods to solve Eilenberger equa- 
tions for the vortex lattice. The main problem in the 
numerical procedure is that the initial conditions for the 
differential equations are unknown. One way to over- 
come that problem is to use special, divergent, gauge 
in which Green's functions are periodic, and solve equa- 
tions in Fourier space (periodic boundary condition).— 
The other method is based on the fact that during the 
integration process Green's functions exponentially grow 
(explode). Fortunately, the exponentially growing un- 
physical solutions can be manipulated to form the phys- 
ical one. This is the essence of the so-called "explosion 
method" iSii^'i^ Here we will use a different approach. It 
is interesting that if one parameterizes quasiclassical Eqs. 
in the form of Riccati's differential equation, then during 
the numerical integration the physical solution is stabi- 
lized regardless of the initial condition. Here we will give 
more details. 

For the s-wave superconductor in presence of nonmag- 
netic impurities Eilenberger equations are 



[w-f u(V + iA)]/ = % + Fg-G/, 



(2) 



u (V - iA)] p = + F*g - Gp. 



(3) 



These are supplemented by the self-consistency equations 
for the gap function ^' and vector-potential A 



^-Ini 2t^ 



UJ>0 



(./>- 



2t 

Vx Vx A = -:^Im^(u5), 



iu>0 



as well as for the impurity potentials 



F 



1 



(/>, 



G 



1 



;.9). 



(4) 



(5) 



(6) 



Born approximation is assumed in treating scattering 
on impurity. For convenience, equations are written in 
following dimensionless units: order parameter is mea- 
sured in units ttTc, length in units Rq — v/{2'kTc), v is 
Fermi velocity, magnetic field in units Hq = <I>o/27ri?Q, 
where $o is flux quantum. Vector-potential is in units 
= <i>o/27ri?o, energy in units i?o = {■kTc)'^NqRI. Scat- 
tering time T is in units l/(27rTc). It can be expressed via 
electron mean free path I: t = ^/0.882^o, with being 
BCS coherence length. Eilenberger parameter k is the 
only material constant that enters the equations 



$0/ (ttTc 



2 ■ 



(7) 



It is related to GL parameter k via — (7C(3)/18)k^ in 
3D case and 



^2 



7C(3) o 

T rv 



(8) 



in 2D case. Here C is Riemmann's zeta function, lu = 
t{2n+l) is Matsubara frequency with integer n, t = T/Tc 
is reduced temperature, u is unit vector directed along 
Fermi velocity. Eilenberger Green's functions /, and 
g are normalized so that g = \/l — ff^. Fermi surface is 
assumed to be isotropic and two-dimensional Average 
over the isotropic cylindrical Fermi surface reduces to 
{■■■) — (I/Stt) / • • • dip, average over polar angle (/?. 

The quantity of our interest, LDOS as a function of po- 
sition r and quasiparticle excitation energy E, is defined 
as 

N{r,E) ^ Na{Re g{r,u,Lu — > S - iE)) , (9) 

where Eilenberger function g describes normal excita- 
tions and (5 is a small number. It is very convenient to 
introduce auxiliary functions a and b through the follow- 
ing transformationiSii^ 



/ = 



2a 



1 + a6' 



2b 



1 + ab 



1 - ab 
1 + ab' 



(10) 
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Equations for auxiliary functions a and b are decoupled 
and have the form of Riccati's differential equation 



uVa = — {uj + G + iuA) a - 



^1 + F 



F*), 
(11) 



+ F* 6^ 

uV6 ={uj + G + iuA) b + y (* + F) . (12) 

Auxiliary functions a{r,u,Lu) and b{r,u,Lu) are not in- 
dependent. Once we solve the Eq. Hll|) . function b can 
be readily calculated: 



b{r,u,uj) = —a*{—r,u,uj). 



(13) 



In the coordinate system (p, ry), where Fermi velocity di- 
rection u coincides with p-axis, 



p ~ X cos (p + y sm <? 
r] — y cos 4> — X sin i 



Eq. (fTT)) reduces to 



iuA) 0.+ 



(14) 



(15) 



A. Iteration procedure 

Iterative procedure for solving Eq. H15|l is the follow- 
ing. We start from some potentials ^'(r), A'{r), F and 
G. It is usual to start from the Abrikosov solution for 
\I'(r), A'{r) = 0, and local values of impurity potentials: 



1 

F = -- 



+ \^{r)\ 



(18) 

After solving the Eq. (|15|) the new values of potentials 
are obtained from the self-consistency eqs. lO, and 
(0, and the new potentials we plug again into the Eq. 
(|15|l and solve it. This iterative procedure is repeated 
until the selfconsistency is achieved. The maximum fre- 
quency Wcut = t{2Ncut + 1) should be chosen so the result 
does not depend on the number of Matsubara frequen- 
cies. On the other hand the number of iteration cycles 
needed to stabilize pair potential increase with the N^ut- 
We followed Klein^^ and choose ujcut = ZOttTc (in physi- 
cal units) as appropriate for various temperatures. This 
gives the number of Matsubara frequencies 



Nr. 



Int 



(19) 



Fortunately it is not necessary to solve Eq. H15|) for all 
Lu. For high frequencies the solution can be well approx- 
imated by: 



where ^' = 'i' + F. Integrating along the direction p from 
p' — Poo to the desired point p' , the physical solution 
a+ is stabilized. Note that integrating in the opposite 
direction, toward decreasing p, one will get solution a_ = 
— l/a_|_. How long integration path p^o should be taken 
depends on lo (is it real or complex) and on impurity 
concentration. 

Vector-potential is written as 



A{r) = h A'{r) 



(16) 



where B is magnetic induction, and A' is periodic with 
V • A' = 0. Therefore, the selfconsistent equation for 
vector-potential can be written as 



Of 

V^A'^^lu,Y,{ug). 



(17) 



It can be easily solved in the Fourier space. 

Auxiliary function a has the same symmetry proper- 
ties as Eilenberger function /, which are described in Ref. 
0. The equilibrium vortex lattice structure is assumed to 
be hexagonal. Therefore, it is sufficient to solve equation 
(|15|1 in the whole vortex lattice cell and only for veloc- 
ity directions < </? < 7r/6. With the help of symmetry 
properties, a(r,M,a;) can be obtained for all velocity di- 
rections. 



1 f 1 



a « — 
2 



{:^ + F), (20) 



where u)' = uj + I/t and 11 = V + iA. For aU n > Ncutl'^ 
we use the equation H2UII . Solution is quasi-periodic. 
Translation by Rnm = nri + mr2 will amount in phase 
factor exp(ix(r, Rnm)) 

a(r + Rnm,v,uj) = a{r,v,uj)e'^^'^'''-\ (21) 

where ri and r2 are primitive vectors of vortex lattice, 
n, m are integers, and 



X = TT 



mx y{n + mcos(3) 



ao 



qq sin /3 



(22) 



The angle between primitive vectors is denoted as /3 
{(3 — 7t/3 in our case). Once the selfconsistent potentials 
^'(r) and A{r) are calculated, the Eilenberger Eqs. are 
solved again but this time for uj = 5 — lE where 5 is small 
number and E is quasiparticle energy. Note that in the 
presence of impurities the Eq. H15(l has to be solved self- 
consistently with respect to impurity potentials F and G. 
As for the choice of 5 one should be very careful. It was 
already noted^fliSi that density of states N{E = 0) is very 
sensitive to the absolute value of 5. Finite 5 has roughly 
the effect of impurities and suppresses the peak in DOS 
at the vortex center. For small values of 5, N{E = 0) is 
spiked at the vortex centers and very fine mesh is needed 
to evaluate average LDOS. We find that S = 0.001 suffice 
for our calculation. 



III. LOCAL DENSITY OF STATES AND 
DIFFERENTIAL CONDUCTANCE 



The physics of the vortex core in the clean Umit is 
very different from the physics of the vortex core in dirty 
superconductors. Properties of the vortex core are gov- 
erned by Andreev bound states in the clean limit, while 
in the dirty limit properties of the core are governed by 
normal electronsiSS To understand the role of impurities 
we briefly explain the formation of bound states. An- 
dreev scattering from the pair potential (order parame- 
ter) inside the vortex core converse electron-like excita- 
tion into hole-like excitation and vice versa. States inside 
the core are coherent superposition of particle and hole 
states. At certain energies the coherent superposition of 
particle and hole states is constructive and the bound 
state is formed. The lowest bound state has the energy 
E = A/kpS,- In the quasiclassical limit kp^ ^ 1, the 
lowest bound state energy is pushed to zero. Zero-energy 
bound state inside the vortex core will manifest as a peak 
in zero-energy LDOS at the vortex center. Scattering on 
impurities will randomize the motion of electron, and the 
coherency is lost. Thus, the impurities will smear out the 
sharp structure of LDOS. To illustrate this we focus on 
the spatial structure of zero-energy LDOS N{r, E = 0). 

In Fig. n spatial variation of LDOS along the line 
connecting two nearest neighbor vortexes is shown. Data 
for a clean superconductor {(,o/i — 0.0), for a relatively 
large mean free path ^o/^ — 0.1 and for impure case 
£,o/£ = 4.0 are presented. Distance between vortexes is 
normalized so that and 1 on the abscissa are position of 
vortexes. To remind the reader again, in the clean limit 
the height and width of the LDOS peak depend on the 
small parameter S, which measures how far we are from 
the pole of the Green function g. In this sense height and 
width of the peak in the clean limit are arbitrary. 

At the vortex core zero energy DOS (ZEDOS) in the 
clean limit highly exceeds the normal state value Nq. 
This was in the beginning at odds with generally ac- 
cepted naive picture of vortex core as being "normal" . 
Analyzing the ZEDOS in the impure case, it is clear that 
coherency is crucial in forming the main peak at the vor- 
tex. In the dirty limit ^o/^ oo ZEDOS within the 
vortex core approaches normal state value iVo, and only 
in this limit one can view vortex core as being "normal" . 
Even a small impurity concentration has a great impact 
on ZEDOS profile. The comparison of the ideal case of 
a clean superconductor Co/^ ~ with rather pure super- 
conductor ^o/^ = 0.1 reveals a change of the vortex core 
size by a factor 2. The change of the vortex core size is 
compensated by the reduction of the peak height, so the 
ZEDOS averaged over vortex lattice cell is approximately 
the same in all cases. 

It is instructive to see how the spatial structure of ZE- 
DOS within the vortex lattice evolves by adding impuri- 
ties. In the clean limit ZEDOS around the single vortex 
is cylindrically symmetric. As soon as vortex lattice is 
formed cylindrically symmetric ZEDOS transforms into 




Nearest neighbour direction 



FIG. 1: Spatial variation of zero energy DOS along the nearest 
neighbor vortex direction. Full line corresponds to the clean 
limit and dashed lines correspond to the superconductors with 
^o/£ — 0.1 and ^o/£ = 4.0. The calculation is performed at 
approximately the same relative field B = 0.1Hc2- 





FIG. 2: ZEDOS within the vortex lattice for superconductors 
with ^o/e = 0.0, ^o/e = 0.1 and Co/^ = 4.0 (in order from left 
to right). Only data points N{E = 0)/No < 1 are presented. 
Small parameter 5 — 0.03 is used for clean limit data to clarify 
the spatial distribution. 



the star-shaped structure within the hexagonal vortex 
lattice,^ This is presented in Fig. [5] The effect of vor- 
tex lattice notwithstanding, the other effects such as the 
anisotropy of the pairing function^^ and the anisotropy of 
the Fermi surface in hexagonal crystal can also contribute 
to the specific star-shaped structure of ZEDOS. By re- 
ducing the mean free path, star-shaped structure gradu- 
ally disappears and is completely absent in the dirty limit 
even at relatively high fields. This indicates that period- 
icity of the order parameter is not the key element to 
explain the structure of N{r, 0) in Fig. |21 Only coherent 
superposition of electron and hole states in the periodic 
vortex lattice can account for the star-shaped ZEDOS. 

In Fig. I2K) LDOS at the vortex center is plot- 
ted as a function of quasiparticle excitation energy 
E (in units ttTc) for the clean case. LDOS oscil- 
lates with energy, the result previously reported in 
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FIG. 3: a) LDOS N{E,r = 0)/A'^o at the vortex center as 
a function of excitation energy E (in units ttTc). b) DC 
a{E)/aN at the vortex center at T = O.lTc- Both data are 
for superconductor in the clean Umit. 



Ref. .24, This phenomenon has the same origin as os- 
cillation of DOS in superconducting-normal proximity 
systems interference of quasiparticles reflected at 
the superconducting-normal barrier. The mixed state 
can be viewed as periodically arranged infinite number 
of "normal" -superconducting boundaries. Here the vor- 
tex cores play the role of normal region in the sense that 
gap drops to zero at the vortex axes. In Fig. Ob) DC at 
T = O.lTc, calculated according to Eq. (Q), is presented. 
At this temperature DC is thermally broadened LDOS, 
but the oscillating pattern is still visible. 

The coherency of quasiparticles is essential both for 
zero-bias peak and oscillation of LDOS with energy at 
the vortex axis. In Fig. QlLDOS at the vortex center as a 
function of energy is plotted for various values of impurity 
concentration. Oscillation amplitude is very sensitive to 
the presence of impurities and is almost lost even in very 
clean samples with S,o/i = 0.1. Proliferating impurity 
concentration will manifest as a flattening of LDOS at 
the vortex center: disappearance of zero-energy peak of 
LDOS, as well as disappearance of deep minima for E < 
^'(B = 0). In the dirty limit ^o/^ ^ oo, LDOS at the 
vortex center is equal to A^o for all quasiparticle energies 
N{r = 0,£:) = NoMii^ 



IV. LOCAL DENSITY OF STATES AND 
SPECIFIC HEAT 

The low energy quasiparticle excitations play the im- 
portant role in the low temperature thermodynamics. 
Specific heat Cs{T) of a superconductor is given by 



T 



dE 



dN{E) 
dT 



In 



2 cosh 



E 

2r 




1 


1 








-3 -1^1 3-3 -1^1 



FIG. 4: LDOS at the vortex center as a function of energy 
plotted for various values of mean free path I. Normalized DC 
o /gn is almost indistinguishable from the normalized LDOS 
N{E,r = Q)/No at this temperature T = O.lTc. 



E , E 
—— tanh — - 

2r 2r 



E^ N{E)AE 




" E ' 


cosh 


2T 





(23) 



One can utilize this expression only if the energy depen- 
dent, spatially averaged, LDOS N{E) is provided. How- 
ever, in the limit T — + the first integral is zero. For 
small T the function to be integrated in the second in- 
tegral is nonzero only in the small vicinity of = 0. 
Therefore we can replace N{E) by N{E = 0) 



Cs 

fim — 2 

T^o T 



E^ N{E = Q)dE _ 2n-'N{E = 0) 





■ E ' 




cosh^ 


2T 









(24) 

In the normal phase Cn/T = 2tt'^No/3 which gives us the 
well known result 



T™ Cn No 



(25) 



If the low energy quasiparticles are localized in the 
vortex cores, which is true for s-wave superconductors at 
least in the limit of very small fields, then N{E = 0) 

I SccW- Here S'ccii = ^o/B is vortex lattice cell area 
and p is the size of the vortex core. If we further as- 
sume that /O^ ~ ^o/Hc2 then we arrive at the following 
scaling relationship N{E = 0) ^ B/Hc2, for s-wave su- 
perconductors. However there is a number of reports of 
nonlinear field dependence of 7s(iJ) in s-wave supercon- 
ductors. One of the offered explanations is that vortex 
core size p itself is field dependent which in turn lead 




to the nonlinear field dependence of zero-energy DOS. 
The shrinking of the vortex core with increasing field is 
detected in NbSe2^ and YBazCugOg.eo— This is fur- 
ther supported by numerical calculations in dirtjii^*^ 
and clean^ limit. Such an explanation brings out an- 
other puzzle. Experimental study on influence of non- 
magnetic impurities on the 7(JJ) in Y(Nii_a;Pta;)2B2C 
and Nbi_a;Taa;Se2 revealed that linear ^{H) is achieved 
only in dirty samples.'^ This result suggest that the vor- 
tex core size in the dirty superconductors is field indepen- 
dent. Numerical calculation by Golubov and Hartman,^'^ 
as well as Sonnier et alpS- shows quite contrary, that even 
in the dirty limit p should shrink with increasing field. 

Here we emphasize the necessity to perform the calcu- 
lation at low temperature in order to analyze the specific 
heat data through the ZEDOS. In Rcf. 8 calculation per- 
formed at r = 0.5Tc revealed that N{E = 0) - C^{B)B, 
where S,{B) is independently calculated vortex core ra- 
dius. At lower temperatures, due to the Kramer-Pesch 
effect, the core radius is smaller and it might have differ- 
ent field dependence. 

The result for the field dependence of ZEDOS in the 
clean limit, for T = O.lTc, is shown in Fig In the 
inset we plot the field dependence of the core radius at 
the same temperature. We define the core radius ^ as 
1/C = {d\'i'{r)\/dr)/\'i'NN\ where I^'atatI is the maximum 
value of the order parameter along the nearest neigh- 
bor direction, and derivative is taken along the same 
direction. Compared to the previously reported result 
at higher temperature T — 0.5Tc, where £,{B) decreases 
with fieldf^ at T = O.lTc, vortex core radius is rather con- 
stant at low fields. As a consequence, zero-energy LDOS 
is also linear function of magnetic induction. 

In the clean limit ZEDOS in between vortexes is neg- 
ligible in fields as large as i? = 0AHc2- In other words. 



the main contributions to ZEDOS is coming from the 
vortex cores. On the other hand, in the dirty limit, 
ZEDOS is not confined to the vortex cores, but it is 
spread throughout the vortex lattice cell. It is large 
even in between vortexes. Thus, the scaling relation 
N{E = 0) ^'^{B)B is of no use in the dirty limit. This 
is the reason why wc do not attempt to correlate vor- 
tex core size £,{B) and field dependence of LDOS in the 
impure case. However, N{E = 0, B) is a linear function 
of magnetic induction at low fields for any impurity con- 
centration: N{E = 0, B) /No = c{t)B/Hc2- Constant of 
proportionally c(t) weakly depends on the electron mean 
free path and saturate s to c 0.8 in the dirty limit. Nu- 
merical calculation of N{E — 0,B)/No as a function of 
mean-free path value is presented in Fig. We note the 
concave curves for dirtier cases. This behaviors coincide 
with the analysis near Hc2 by Kitai^ 

In Fig. [7| is shown the field dependence of the core 
radius as calculated from the pair potential profile ^'(r). 
For a fixed relative field B /Hc2 core radius ^ is a nomono- 
tonic function of mean free path, first sharply increases 
and then slowly decreases with increasing of ratio In 
the dirty limit vortex core shrinks with increasing field, 
which is consistent with the previous calculations^^sS^ 
in sharp contrast with vortex core enlargement with in- 
creasing field in the clean limit. 

The experimental data, however, revealed that con- 
stant c ~ 1 in the dirty limit. It also shows that scaling 
N{E = 0,B)/No = c{t)B/Hc2 for all field values, which 
is a remarkable feature that still lacks the explanation. 
Worth is mentioning that in Ref . specific heat is a non- 
linear function of field in samples Y(Nii_2;Pt j.)2B2C for 
all < a; < 1. In these materials, we need to consider 
also the effect of gap anisotropy. 
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that coherency is crucial in forming the spatial structure 
of LDOS. As soon as impurities are introduced into the 
superconductor, scattered electrons loose the information 
on their initial state, the coherency is lost and sharp 
LDOS structure is flattened. It is calculated how DC 
spectra evolve with electron mean-free path. Although 
the impurities have a great impact on LDOS, spatially 
averaged LDOS shows weak dependence on relative field 
B/Hc2- We hope that present calculation can be helpful 
to roughly estimate the electron mean free path through 
the LDOS measurement. 



FIG. 7: Field dependence of vortex core size for various mean 
free path. 
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